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Abstract 

We calculate the differential cross section and single-polarization observables 
for the reaction jp ^ ^p within perturbative QCD, treating the proton as 
a quark-diquark system. The phenomenological couplings of gauge bosons to 
(spatially extended) diquarks and the quark-diquark distribution amplitude 
of the proton are adopted from previous investigations of baryon form factors 
and two-photon processes. Going beyond leading order, we take into account 
hadron-mass effects by means of a systematic expansion in the small parame- 
ter (hadron mass/ photon energy). With the <&-meson distribution amplitude 
taken from the literature our predictions for the differential cross section at 
|t| ^ 4GeV^ seem to provide a reasonable extrapolation of the low-t data and 
are also comparable in magnitude with the results of a two-gluon exchange 
model in which the gluons are considered as a remnant of the pomeron. For 
momentum transfers of a few GeV hadron-mass effects appear still to be size- 
able. 
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I. INTRODUCTION 

In the course of the last few years exclusive photo- and electroproduction of neutral 
vector mesons (p, cu, $, J/^P) at intermediate and high energies has attained increased 
theoretical interest, which has been stimulated by corresponding experimental efforts at 
CEBAF [|T| and HERA [Q. The aim of the HERA experiments is a more fundamental 
understanding of Pomeron phenomenology in terms of QCD. Vector-meson dominance in 
combination with Pomeron exchange proves to be very successful in describing diffractive 
photo- and electroproduction of neutral vector mesons, provided the energy is the only large 
scale in the process f^. In their approach Donnachie and Landshoff |Q assume that the 
pomeron couples directly to on-shell quarks in the vector meson and the proton, respectively. 
A QCD-inspired version of Pomeron exchange, in which the Pomeron is replaced by two 
non- perturbative (abelian) gluons, goes back to the same authors and is usually termed 
a "soft Pomeron". In the limit t ^ the soft-pomeron exchange becomes very similar 
to the usual Pomeron exchange. This can be used to interpret the form factor at the 
quark- Pomeron vertex [Q. However, if in addition to the energy another scale, like the 
photon virtuality or the mass of the vector meson, becomes large perturbative QCD effects 
tend to enter the game. For large photon virtuality (and small values of Mandelstam t) 
perturbative QCD becomes applicable due to a factorization theorem which has been proved 
in Ref. |^. The three building blocks of the factorization formula are the hard photon- 
parton scattering amplitude, the distribution amplitude of the vector meson and finally a 
skewed gluon distribution of the proton. As a link between hard inclusive and exclusive 
reactions which generalizes the ordinary parton distributions, skewed parton distributions 
have recently attracted a great amount of interest (see, e.g., Ref. and references therein). 
Electroproduction with highly virtual photons thus offers a possibility to study such skewed 
parton distributions. Within this perturbative formalism for vector-meson electroproduction 
with highly virtual photons the two-gluon aspect of the pomeron emerges automatically in 
a form which is related to the skewed gluon distribution in the proton. 



One of the aims of photoproduction experiments at CEBAF is, on the other hand, to 
understand the transition from non-perturbative to perturbative production mechanisms. 
In the case of photoproduction the transverse momentum transfer p\ ^ tu/ s has to be 
large in order that perturbative QCD becomes apphcable. The pertinent factorization for- 
mula, which in general holds for arbitrary exclusive hadronic processes at large transverse 
momenta, is the result of a leading-order perturbative analysis. It represents a hadronic am- 
plitude as a convolution integral of a hard-scattering amplitude with distribution amplitudes 
(DAs) 1^ . The process-dependent hard-scattering amplitude is perturbatively calculable and 
represents the scattering of hadronic constituents in coUinear approximation. The process- 
independent DAs contain the non-perturbative bound-state dynamics of the hadronic con- 
stituents. There is no doubt that this hard-scattering approach (HSA) gives the correct 
description of exclusive hadronic processes in the limit of asymptotically large momentum 
transfer. But there is increasing evidence that the pure perturbative contribution is not the 
dominant one at experimentally accessible momentum transfers. A reasonable reproduc- 
tion of the existing exclusive large-t data is only achieved with very endpoint- concentrated 
hadron DAs, like those proposed by Chernyak and Zhitnitsky on the basis of QCD sum 
rules [0. Such DAs, however, occur to be problematic for a perturbative calculation since 



they enhance contributions from the soft end-point regions in the convolution integral [[10 
Some improvement with regard to the consistency of the perturbative calculation can be 
achieved by applying a modified factorization scheme in which transverse-momentum effects 
and Sudakov suppressions are included. But this leads to a substantial reduction of the per- 
turbative contribution [|lT|. For the pion an end-point concentrated DA meanwhile seems 



even to be ruled out by the recent CLEO data [|T^] for the 7r7-transition form factor, since 
it provides a much too large perturbative result. On the other hand, it has been shown 



recently for electromagnetic nucleon form factors |[T^ and for wide-angle Compton scatter- 
ing ||7|,p!^ off protons that the existing data can be saturated by soft overlap contributions. 
Furthermore, these contributions can be modeled by a nucleon light- cone wave function 
which gives rise to a distribution amplitude close to the asymptotic one. This, however. 



implies that the perturbative contribution becomes small. 

The approach which we are going to apply for hard exclusive photoproduction of $ 
mesons is still based on the HSA but includes also a certain amount of soft physics which is 
modeled by means of diquarks. Diquarks may be considered as an effective way to account for 
quark-quark correlations in baryon wave functions. Part of the soft overlap contributions are 
then absorbed into phenomenologically parameterized diquark form factors which occur at 
the gauge-boson diquark vertices. This HSA-based diquark model has already been apphed 
successfully to other photon-induced hadronic reactions like magnetic and electric baryon 
form factors in the space- |T^ and time-like region |T^ , real and virtual Compton scattering 



17,18], two-photon annihilation into baryon-antibaryon |1|,1S], and photoproduction of K 



mesons 



T9[| . Further applications of the diquark model include the charmonium decay rjc 



pp [jT6[ and the calculation of Landshoff contributions in elastic proton-proton scattering 
pO| . The purpose of the present investigation is to extend this foregoing work to 7 p -^ $ p. 
The reason why we concentrate on the $ channel is that it is cleaner and simpler than 
photoproduction of p or cj mesons. Within the HSA the valence Fock state of the $, i.e. 
the s-s state, can only be produced via the exchange of two gluons between the $ and 
the proton. Quark exchange could only happen via the knockout of an s-s pair from the 
strange-quark sea in the nucleon which is, however, power suppressed within the HSA. 

In the following section we start with a short outline of the diquark model, show 
how Feynman diagrams are properly grouped according to their propagator singularities 
(Sec. [II B|) , sketch how we deal with the propagator singularities numerically (Sec. |II C| ) , 
and present our treatment of hadron-mass effects (Sec. [II D| ). Section |ITI| summarizes how 
photoproduction of vector mesons is described in terms of helicity amplitudes. Furthermore, 
the fixed-angle scaling behavior of these helicity amplitudes is discussed and related to the 
treatment of mass effects. Analytical expressions for the hard part of the hadron-helicity- 
conserving amplitudes are also given in this section. Our predictions for photoproduction 
observables along with predictions from Pomeron and soft-Pomeron exchange are presented 



in Sec. |V| At the end of this section we comment on the prospects of extending our calcu- 



lation to photoproduction of J/\l/ mesons. Finally, a summary and concluding remarks can 
be found in Sec. |V|. 

II. FORMALISM 

A. The hard-scattering formalism with diquarks 

As already mentioned before, within the HSA an exclusive scattering amplitude at large 
momentum transfer is expressed as a convolution integral of a hard-scattering amplitude 
with hadron distribution amplitudes [Q. For the photoproduction reaction 7p — > ^p, 
which we are interested in, this integral takes on the particular form 

M{x}{s,i)=J dxidyidzicj)^ {zi,Pi.)(t)^\yi,p±)T{x}{xi,yi, Zi; s,t)(f)P{xi,p±_) . (2.1) 

The distribution amphtudes 0^ are probability amplitudes for finding the valence Fock state 
in the hadron H with the constituents carrying certain fractions of the momentum of their 
parent hadron and being collinear up to a factorization scale p_j_. In our model the valence 
Fock state of an ordinary baryon is assumed to consist of a quark (g) and a diquark (D). We 
fix our notation in such a way that the momentum fraction appearing in the argument of 0^ 
is carried by the quark - with the momentum fraction of the other constituent (either diquark 
or antiquark) it sums up to 1 (cf. Fig. |l|). For our actual calculations the (logarithmic) 
pL dependence of the DAs is neglected throughout since it is of minor importance in the 
restricted energy and momentum-transfer range we are interested in. The hard scattering 
amplitude T{a} is calculated perturbatively. Intrinsic transverse momenta of the constituents 
are thereby neglected, i.e., one assumes that the momenta of the constituents are collinear 
to those of their parent hadron. T^x} consists in our particular case of all possible tree 
diagrams contributing to the elementary scattering process 'y q D -^ QQqD. Eight of the, 
altogether, sixteen diagrams which contribute to photoproduction of (heavy) quarkonia are 
depicted in Fig. |^. The subscript {A} of T represents a set of possible photon, proton and 
vector-meson helicities. For our purposes it appears to be more convenient to express the 
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analytical results in terms of massless Mandelstam variables s, i, and u than in terms of the 
usual (massive) ones. A detailed account of our treatment of hadron-mass effects is given 
in Sec. [II D| below. 



The diquark model, as applied in Refs. [0- ||T9[, comprises scalar (S) as well as axial 



vector {V) diquarks. The dynamics of S and V diquarks is determined by their coupling 
to photons and gluons. The expressions for the photon-diquark and gluon-diquark vertices 
correspond (almost) to the most general form for the (parity and time-reversal invariant) 
coupling of a spin-1 gauge boson to a spin-0 or a spin-1 particle, respectively |T^. Thereby V 
diquarks are allowed to possess an anomalous (chro mo) magnetic moment Ky- In applications 
of the model Feynman graphs are calculated first with the corresponding Feynman rules for 
point-like diquarks. As a second step, Feynman graphs with {n — 2) gauge bosons coupling 
to the diquark are combined in a group which we call "n-point contribution". Finally the 
composite nature of diquarks is taken into account by multiplying each n-point contribution 
with phenomenological vertex functions, the diquark form factors. The particular choice 

F^iQ') = Ss ,S'r... ^ (2-2) 






n'\Q')=5y\7i^^] , (2-3) 



for 3-point contributions and 

F(") = a5Ff(Q^), (2.4) 



)2 \ {"~3) 

Fr = avF^'\Q'')\T,P^] , (2.5) 



-W _ . I?(3)^n2^ { _Qv_ 

.Ql + Q' 



for n-point contributions (n > 4) ensures that in the limit Q^ ^ oo the scaling be- 
havior of the diquark model goes over in that of the pure quark HSP. The factor 
^s(v) = ci:s{Q^)/c(s{Q%v)) i^s{v) = 1 for Q"^ < Q'^s(v)) Provides the correct powers of 
as{Q^) for asymptotically large Q^. For the running coupling a^ the one-loop result 
as = 127r/251n((5^/AQ(^^) is used with Aqcd = 200MeV. In addition as is restricted 
to be smaller than 0.5. as and ay are strength parameters which allow for the possibility of 
diquark excitation and break-up in intermediate states where diquarks can be far off-shell. 



Having sketched how diquarks are treated perturbatively, a few words about the choice 
of the DAs (which incorporate the bound-state dynamics of the hadronic constituents) and 
the spin-flavor wave functions of the proton and the $ are in order. If one assumes zero 
relative orbital angular momentum between quark and diquark and takes advantage of the 
collinear configuration {pg = XiPb and po = X2PB = (1 ~Xi)pb) the valence Fock-state wave 
function of an octet baryon B may be written as 

^^(pb; A) = /#0f (xi) xs uiPB, A) + /# 0^(xi) xv 4 (7" + — ) iMpb, A) . (2.6) 



The two terms in Eq. ( p.6|) represent states consisting of a quark and either a scalar or a 
vector diquark. The pleasant feature of the covariant wave-function representation, Eq. ( p. 61) , 
is that it contains, besides Xi and a (the Lorentz index of the vector-diquark polarization 
vector), only baryonic quantities (momentum pb, helicity A, baryon mass tti^). Assuming 
an SU(6)-like spin-flavor dependence, the flavor functions Xd for a proton take on the form 

X^s = uSiu,d], Xv= [uV{u,d}-V^dV{u,u}]/\^- (2.7) 

In Eq. ( |2.6| ) the color part of the quark-diquark wave function has been neglected. Since 
the quark-diquark state has to be a color singlet the diquark has to be in a color antitriplet 
state (like an antiquark) Thus the quark-diquark color wave function is simply ipqif^ = 

(1/v^) ELi Sa-a. 

Like the baryon wave function, Eq. ( p.6|) , also the Q-Q wave function of a vector meson 
M can be represented in a covariant way 

^""(Pm; A) = -f' 0^(^i, A) x^' ^{i>M + ^m) /(A) , (2.8) 

with e(A) denoting the polarization vector of the vector meson. The flavor function of the 
$ is X* = ss. 



In Refs. [0- [|T^ a quark-diquark DA of the form {c\ = C2 = for S diquarks) 

(2.9) 



0^(a;i) = A^-"xi(l - xi)^(l + CiXi + c^xl) exp 



\Xi (1-xi, 



turned out to be quite appropriate for octet baryons B. The origin of the phenomenological 
baryon DA, Eq. ( |2.9| ), is a non-relativistic harmonic-oscillator wave function (cf. Ref. ^^^)- 
Therefore the masses appearing in the exponential have to be considered as constituent 
masses (330 MeV for light quarks, 580 MeV for light diquarks, strange quarks are 150 
MeV heavier than light quarks). The oscillator parameter 6^ = 0.248 GeV~^ is chosen 
in such a way that the full wave function gives rise to a value of 600 MeV for the mean 
intrinsic transverse momentum of quarks inside a nucleon. The exponential does not only 
provide a flavor dependence, but also suppresses the end-point regions xi, yi — > 0, 1 in the 
convolution integral, Eq. (|2.1| ). One thus avoids to pick up sizable contributions from the 
soft end-point regions and results become less dependent on details of the baryon DA. In 
the actual data fitting the exponential plays only a minor role. Therefore, the constituent 
masses and the oscillator parameter have not been considered as free parameters but have 
been fixed in advance. We stress that the constituent masses do not appear in the hard 
scattering amplitudes (cf. Sec. [11 D| below). 

A model for the $ meson DA which fulfills QCD sum-rule constraints has been proposed 
by Benayoun and Chernyak ||22|. Since we also want to investigate to which extent our 
calculations for $ production may be generalized to J/\l' production we modify this DA by 
attaching an exponential factor which provides a flavor dependence. The effect of this factor 
is negligible for $s, but it becomes crucial for J/\l's. The form of the exponential factor 



is inspired by a simple relativistic treatment of heavy mesons [^ and provides a flavor 
dependence in accordance with heavy-quark effective theory. For longitudinally polarized 
$s and J/\E's the DA reads (M = $, J/\E') 

/^(zi, A$ = 0) = N^i^o^iil - zi) exp [-Pml,{z, - 0.5)^] . (2.10) 

The DA for transversally polarized $s and J/'^s contains an additional polynomial in zi 
which enhances the maximum at zi = 0.5 and makes it narrower (cf. Ref. p2|) 



r'{z,, X^ = ±1) = N;^i=^,ziil - z,) exp [-b'mi,{zi - 0.5) 

X (5225.2 zl{l - zif + 0.39 zi{l - ^i)) . (2.11) 



For the oscillator parameter b we take a value of 0.97 GeV ^ in accordance with estimates 



for the radius of the J/'^ meson [|J] 



The usual normalization condition, /q dx(j)^{x) = 1, fixes the constants A^'^ and A^^ 



in Eqs. (g^), (gAOD , and (|TT|). The quantity /*^ showing up in Eq. (U) is related to 



the experimentally determinable decay constant of the meson M. In case of the $ we get 
from the radiative decay width of the $ -^ e+e^ decay /* = /j^cay/^v^ = 46.95 MeV. 
The analogous constants fg and fy for the q-D DAs of baryon B are free parameters of 
the diquark model. They are, in principle, determined by the probability to find the q-D 
state {D = S, V) in the baryon B and by the transverse-momentum dependence of the 
corresponding wave function. 

Our present study is performed with the set of diquark-model parameters found in 



Ref. |]T5[ by fitting elastic electron-nucleon scattering data. The numerical values of the 
parameters are: 

/_5 = 73.85 MeV, Q| = 3.22 GeV^ as = 0.15, 
fv = 127.7 MeV, Q^ = 1.50 GeV^, ay = 0.05, 
Kv = 1.39, ci = 5.8, C2 = -12.5 . (2.12) 



For further details of the diquark model we also refer to the paper of Jakob et al. [jT5|. 



B. The general structure of the 
hard-scattering amplitude 

The introduction of diquark form factors already requires a decomposition of the hard- 
scattering amplitude Tj^} into 3- and 4-point contributions 

T{x}{xi,yi, zi, s, i) = esfl-^f\xu VuZi, s, i) 

+ esfl_'^}^\xi,yi,zi;s,i) . (2.13) 
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Cq and e5j^ ^j are the electric charges of the quark q and the S diquark S[u,d] in units of 
Co, respectively. The superscripts indicate whether the photon couples to the heavy quark 
Q in the meson, or the light quark q in the proton with the diquark D = S,V acting as 
spectator (3-point contributions), or whether the photon couples to the S* or l^ diquark in 
the proton (4-point contributions). Note that the only V- diquark contribution is Tr^y . 
The reason is that the charge-flavor factor in front of T|^| vanishes and Tr^J represents 
already a mass correction of second order, which is neglected in our calculation. For the 
numerical treatment of the convolution integral, Eq. ( p.l|) , it is advantageous to further split 
the n-point contributions into two parts which differ by their propagator singularities: 



Tll}^\xi,yi,zi;s,i) 



AQ,D) 



iQ,D)i 



/{A} (a;i,2/i,2;i;s,t) g\^y '{xi,yi, zi; s,t) 



+ 



(gi + le) {qi + ie') (gf + le) (g| + ze' 



cil^D), 



Sq,D), 



/{A} \^i^yi^zi]s,t) 5fj^'| '{xi,yi,zi]s,t) 



T, 



P) 



{x}{.xuyi,zi]s,t) 



)l + ie){ql + ie') 



{9i+ie)iqi + ie') 



^§\ 



^gh 



-X2y2i) 



-X2y2i) 



f{x}{xuyi, zi; s, t) g\x}{xuyi, z^; s, t) 



Fff{-X2y2i). (2.14) 



{al + ^e) (gi + «e') {9I + ^e) {ql + ^e') 
The functions / and g correspond to gauge-invariant subgroups of the 16 Feynman diagrams 
which contribute to the hard-scattering amplitude T{a} (cf. Fig. ^. The diagrams entering 
the gs are related to those entering the /s by interchange of the two gluons. The g~^ and 
g^"^ denote just those quark and gluon propagators which can go on- shell when integrating 
over Xi, yi, and Zi. Explicitly, the propagator denominators read: 



ql = y2Z2S + X2y2i + X2Z2U , g|= y2ZlS + X2l/2^ + X2ZiU , 

ql = yiZ2S + xiyii + X1Z2U , g|= yiZis + xiyj + xiZiu , 



9i = xiy2S + X2yiu ; 



9i=X2yis + xiy2U. 



(2.15) 



As already indicated in Eq. ( |2.14|) propagator singularities are treated by means of the usual 
ie prescription. The resulting imaginary part in the amplitudes is considered as a non-trivial 



prediction of perturbative QCD which is unaffected by long-distance effects [^. As we will 
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see below, it gives rise to non- vanishing predictions for polarization observables which require 
the flip of a hadronic helicity. For analytical expressions of the functions / and g and their 
symmetry properties we refer to Sec. |ITI| . 

C. Numerical treatment of propagator singularities 

The numerical difficulties in performing the convolution integral, Eq. (|2.1|) , are mainly 
caused by the occurrence of propagator singularities in the range of the (xi,?/i,Zi)- 
integration. These give rise to principal value integrals 

1 /1\ . . ..,2 



k'^±ie \k 



pi- ±t7rS{k'). (2.16) 



As explained in some detail in Ref. [|T9| such integrals can still be treated by means of a rather 



fast and stable fixed-point Gaussian quadrature after carefully separating the principal value 
parts and doing the corresponding integrations analytically. Here we only want to give a 
short account of our integration procedure. We first perform the xi-integration for which 
the various contributions to the hard scattering amplitude (cf. Eq. ( |2.14| )) lead to integrals 
of the general form 

In order to simplify notations we have neglected helicity labels and the dependence on the 
Mandelstam variables s and t. Furthermore, the distribution amplitudes have been absorbed 
into the function h{xi,yi,zi). By applying a partial fractioning, exploiting Eq. ( p.l6|) , and 
separating the principal value part Eq. (p.l7|) may be rewritten to give 



^(>'^)_^(^^)+^e[\^x^J Jo kf 

dkj 






dxi 
— same (fci <-> ^2) \ ■ (2-18) 
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Xi = Xi (yi,Zi) and Xi = Xi ^\yi, zi) are the zeroes of the propagator denominators 
kf and fc|, respectively, and e stands for either e or e', depending for which of the two 
the hmit e*^'^ — ^ has been taken. As fong as the two zeroes Xi and Xi ^ (which still 
depend on yi and zl) do not coincide in the (|/i, 2i)-integration region - this is the case 
for the contributions of f^'^'^\ g^'^'^\ f^^\ and g^^^ - the (xi,?/i, 2;i)-integration can be 
performed very efficiently with the help of Eq. ( p^.l8| ). After analytic integration of the two 
principal- value integrals in Eq. ( |2.18| ) (which is a simple matter) all the integrands of the 
remaining integrations are regular functions of the integration variables and can be treated 
by fixed-point Gaussian quadrature. 

For /W'^) and g^^'^\ however, the situation becomes a little bit more complicated. After 
having performed the Xi-integration in Eq. (p.l8|) one observes that the two propagator 



singularities xi and x\ ^ coincide along a trajectory in the (yi, 2;i)-integration region. 
This leads to a quadratic zero yj ^' = yj ^' '^ (zi) of {xi — Xi ). A closer inspection 
of Eq. (|2.18|) reveals that this quadratic zero in the denominator of l'^^^'^^\yi^ Zi) does 



not cause any further problem when performing the (yi, 2;i)-integration for the real part of 
/(fci,'c2)^y_i^^ 2:1). The reason is that it is completely compensated by a corresponding quadratic 



(fci) 



zero in the numerator of 'StI^^^'^'^\yi^Zi). In ^l'^^^'^'^\yi^zi)^ however, the zero of {xi'^' — 
Xi ^ ) is only partly compensated by the numerator and one still encounters a single pole. 
At the pole '^l''^^'^^\yi, zi) becomes proportional to {yi — y| ^' )/[{yi — yi ^' )^ + ie] = 
{1/(1/1 - yf^'^''^ + ie) + 1/(1/1 - i/f^'^''^ - ie)}/2. As a consequence J^ dyiQI^''^'''^\yi, zi) 
becomes a pure principal- value integral (cf. Eq. (|2.16|) ) 

f'dy^^I^''^''^^\y,,z,) = p f'dy^^I^''^''^\y,,z,) (2.19) 

Jo JO 

which can be treated numerically analogous to the principal- value integral in Xi. One only 
has to separate the principal-value part and do the corresponding integration analytically. 
For this purpose the analytical expression of the residue of QI^''^'^^\yi, zi) at the position 
of the pole yj ^' ^ (zi) is needed. The remaining integrations are then again amenable to 
Gaussian quadrature. 
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Thus it is possible, by carefully separating the singular contributions, exploiting S func- 
tions, rewriting principal-value integrals as ordinary integrals plus analytically solvable 
principal- value integrals, to do all the numerical integrations by means of fixed-point Gaus- 



sian quadrature. Whereas the authors of Ref. [^, who have used a different integration 
procedure, were only able to obtain reliable results for the imaginary part of the convolution 
integral, Eq. (|2.1| ), our numerical results for the imaginary and the real part occur to be 
very stable. Going from an xi x yi x z\ grid of 24 x 24 x 20 points to a 48 x 48 x 32 grid 
changes the results for the amplitudes (over the whole angular range) by less than 0.1%. 
Our numerical calculations were performed on an AlphaServer 1000 4/266. For the larger 
grid size the evaluation of da / dt{'-)p — > ^p) (taking into account the 8 hadronic non- and 
single-flip amplitudes) took less than 7 seconds per energy point and angle. 

D. Treatment of mass effects 

Our calculation of the hard-scattering amplitude involves an expansion in powers of 
{rriH/y/I) which is performed at fixed scattering angle. We keep only the leading order and 
next-to-leading order terms in this expansion. Hadron masses, however, are fully taken into 
account in flux and phase-space factors. The reasons why we want to include hadron-mass 
effects in the hard scattering amplitude are twofold. On the one hand, we want to apply our 
model already at momentum transfers of only a few GeV, where hadron masses may still 
play a role. On the other hand, we also want to make predictions for spin observables which 
require the flip of hadronic helicities. Within usual perturbative QCD such observables 
would vanish since a spin-1 gauge boson which couples to a nearly massless (current) quark 
cannot flip the helicity of the quark. In the original diquark model |1^ violations of hadronic 



helicity conservation are generated by massive V- diquarks on the account of introducing a 
mass parameter for the V-diquark. In this paper we adopt another strategy which allows 
for a more consistent treatment of mass effects without introducing new mass parameters 
for the hadronic constituents. Like the authors of Ref. ET^ we assume for every hadronic 
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constituent that its 4-nionientuni is proportional to the 4-nionientum of its parent hadron. 
This requires to assign to every constituent of the hadron H an effective mass xtuh, where x 
is the fraction of the four- momentum of the hadron H carried by the constituent. Keeping 
in mind that the momentum fractions are weighted by the hadron DAs in the convolution 
integral, Eq. (^.11), this means that quarks and diquarks acquire (on the average) a mass 
which is rather the mass of a constituent than a current particle. This may be interpreted 
in such a way that at intermediate momentum transfers a quark or diquark surrounded by 
a cloud of qq pairs and gluons is not resolved but rather acts like a single effective particle 
with a corresponding mass. 

This prescription for the constituent masses applies to the (on-shell) particles which cor- 
respond to the external legs of the Feynman diagrams. It has already been tacitly exploited 
when casting the hadronic spin wave functions into covariant form (cf. Eqs. ( |2.6| ) and (|2.8|) ). 
Taking the quark-diquark DA of Eq. ( p.9|) the average values of the effective masses of quarks 
and diquarks in a proton become 



i'mf) = {xi)mp ^ -TUp , 



1 

-i 

3 



K^fif) = ((1 - x,))mp ^ 'Krip . (2.20) 

Analogously the effective mass of the heavy (anti) quark in the meson becomes on the average 
half of the meson mass. The choice of the masses which occur in the propagators of internal 
particles is not that obvious, since, e.g., a quark in the incoming proton has mass X\mp 
whereas the same quark in the outgoing proton has mass yinip. So what should we take for 
the mass of this quark between two vertices? In order to explain our recipe we, e.g., consider 
the momentum of the heavy quark between the two gluon vertices in the third diagram of 
Fig. |. It is 

q5 = zik + yipf -xipi. (2.21) 

with Pi, pf and k representing the 4- momenta of the hadrons as indicated in Fig. |I|. The 
quark mass in the propagator related to q^ is now obtained by replacing the 4-momenta in 
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Eq. ( |2.21|) by the masses of the corresponding particles, i.e. 

m§'^ = ziniM + (2/1 - xi)mp . (2.22) 

Taking again into account that the momentum fractions are weighted by the hadron DAs, 
Eq. ( |2.9| ) and Eq. ( |2.10| ), the average heavy quark mass in the propagator becomes {ttiq ) ^ 
0.5 mj\/. Extending this recipe to all the quark and diquark propagators and using the 
covariant wave-function representations, Eqs. (|2.6|) and (p.8|) , for the external particles, leads 
to analytical expressions for the Feynman diagrams which contain (apart of the momentum 
fractions Xj, yi, and Zi) only hadronic momenta, masses, spinors, polarization vectors, and 
helicities. The occurring traces and spinor products can be evaluated in the usual way. The 
final step is then to reexpress the (massive) Mandelstam variables s, t, and u in terms of 
massless ones 

s = s(l + 0(m^/s)) with s = 4g^ , 

t = i{l + 0{mjj/s)) with t = -2g2(l - cos(^cm)) , 

M = M(l + 0(m^/s)) with u = -2q^{l + cos{9cm)), 

(2.23) 

and to make an expansion in terms of the (small) parameters {mu/VS), keeping the scat- 
tering angle 6cm fixed. In order to facilitate the treatment of propagator singularities we 
retain only the leading and next-to-leading order term in this expansion. Taking again as 
an example the propagator related to the momentum q^ (cf. Eq. ( |2.21| )) the propagator 
denominator becomes 

(?5 - ^q""^ ) = yi^i(s - {mp + muf) + XiVit + XiZi{u - {nip - mMf) 

= yizis{l + 0{mjj/s)) + Xiyii{l + 0{mjj/s)) + xiZiu{l + C(m^/s)) 
^ql (2.24) 

hence the quark propagator is given by 
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I.e., by neglecting mass corrections of 0{m\j/s) we end up with propagator denominators 
which still appear like those of massless particles (cf. Eq. ( p.l5| )). Apart from the simplifi- 
cations in the treatment of propagator singularities and in the analytical expressions for the 
hard part of the amplitudes our approximate treatment of mass effects has a few other nice 
features. By including mass corrections of C(mj7/-\/I) we obtain non-vanishing results for 
hadron-helicity-flip amplitudes and thus for polarization observables which require the flip of 
a hadronic helicity (see the discussion in the next section). Like the leading-order hadron- 
helicity-conserving amplitudes also the hadron-helicity-fiip amplitudes are U{1) gauge in- 
variant with respect to the photon and SU{'i) gauge invariant with respect to the gluon. 
Finally, crossing relations for the hadronic helicity amplitudes are fulfilled up to 0{mH/y/S), 
as has been shown for Compton scattering off baryons and its crossed process 77 -^ BB |]18| . 



III. HELICITY AMPLITUDES 
AND OBSERVABLES 

For exclusive photoproduction of vector mesons 7P — > M p one finds, altogether, 24 
complex helicity amplitudes. By virtue of parity invariance only 12 of these helicity ampli- 



tudes are independent. Following the notation of Ref. ||28| we denote them by 



-f^l.A* = ^A<[>,Ay=+l/2,A^=l,Ai=-l/2 ; 

-f^2,A<i, = ^A$,Ay=+l/2,A^=l,Ai=+l/2 , 

-f^3,A<s, = ^A$,A/=-l/2,A^=l,Ai=-l/2 ; 

Ha^X^ = ^A4>,Aj=-l/2,A^=l,Ai=+l/2 , (3.1) 

with A$ = 0, ±1. For our normalization of the helicity amplitudes the unpolarized differential 
cross section takes on the form 

"^^ ' E hH.M?. (3.2) 



dt 327r(s-m2)2^^^Q_^^.^^ 
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As examples for single spin observables, for which our model provides nontrivial predictions, 
we will consider the beam and target asymmetry. The corresponding expressions are |^ 

„ da dan da \ 



^ dt dt dt 



167r(s — rrip)^ 



^ (-^^4,1-^1,-1 ^ HIqHi^q + Hl^^Hii — H^^H2-i + H^QH2fi — Hl_^H2,ij (3.3) 



for the photon and 

^da 1 



^ dt 167r(s - m2)2 

X 0= [Hl_^H^^^i + HIqH^^q + Hl^^H^^i + H2_iHi_i + H^qHiq + H^^^Hi^ij (3.4) 

for the proton, respectively. 

Within the hard scattering approach the energy dependence of the helicity amplitudes 
at fixed cm angle and large s is roughly 

H2,o,Hs,o (xs-''/\ 

Hlfi , -^4,0 , -f^2,l , Hs^i, H2-I , -^3,-1 OC S~ , 

depending on whether the helicity of the hadronic constituents is conserved or flipped by one 
or two units. This scaling behavior is modified by logarithms due to the running coupling 
as and, eventually, the evolution of the hadron distribution amplitudes. Further deviations 
are due to the diquark form factors. For small momentum transfers diquarks appear nearly 
point- like and the decay behavior of the helicity amplitudes is weakened by one power of 
s. The energy dependence in Eqs. (|3.5|) is only approached at large enough momentum 
transfer, where the diquark form factors become fully operational. 

For the helicity amplitudes our treatment of mass effects (cf. Sec. [1IE)| ) leads to an ex- 
pansion in powers of {rriH / \/§y . As compared to the hadron-helicity-conserving amplitudes 
7^2,0 and ^^3,0, the leading-order term of the single-flip amplitudes Hiq, H^q, H2^±i, and 



iJs.ii is suppressed by 0{mH/\s). The leading-order term of the double-flip amplitudes 
Hi±i and i/4,±i is even suppressed by {niH/ys)'^ and thus has the same order as the first 
(non- vanishing) mass-correction term of the hadron-helicity-conserving amplitudes. To sim- 
plify matters we only consider non- and single-flip amplitudes to leading-order and neglect 
double-flip amplitudes at all. 

In the following we quote the analytical results for the functions / and g which have been 
obtained with the help of the symbolic computer program "Mathematica" , in particular with 
the two program packages "FEYNARTS" fH] and "FEYNCALC" @ . In case of the hadron- 
helicity conserving amplitude -^3,0, i-e. for the helicity combination A<j, = 0, A/ = —1/2, 



A. 



+1, Aj = —1/2, the functions / read: 



/ 



(Q,5) 



^^l AQ,V) 



{o,-|;+i,-|} ^i '{o.-^;+i'-i} 



-2V- 



u 



Z2-Z1) ziZ2V^{xi{l + X2)~yi{l + y2)) 



f 



i<l,S) 

{o,-i;+i,-^} 



XlX2yiy2ZlZ2t^ 

+iu{x2 + y2){xi - zi){yiZ2 + y2Zi) - y2zlP{^ + ^2) 

+xiy2p 4:zl{yi - zi) + yi{x2 + ^2) + zi{2xiZi - 3|/i) } 

1 



^jK Aq,V) 
1 



-4V- 



X2yit 



{X2Z1U + y2Z2s] 



|u, 2,+i, 2 1 xiyit 



{xiziu + yiZ2s} 



(3.6a) 
(3.6b) 

(3.6c) 



with 
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{o,-|;+i-4} ^{o,-i;+i,-i} 



V = 327i^^/^a1{iu/s)\/~iCF , (3.7) 

where Cp = 2/3-\/3 denotes the color factor and a the fine-structure constant, respectively. 
The functions g are related to the /s in a simple way: 

^^^'"^ -/iri.,. n (-1--2), (3.8a) 

(3.8b) 

si;;_.,,_.} = -/r„:-4.«,-i} (^. - ^j ■ (^.s^) 

Vector 4-point contributions f^^^ and (7*-^'' are, in general (for all helicity combinations), 
suppressed by 0(m^/s) or even stronger. Relations analogous to Eqs. (|3.8|) are also valid 






K5) 
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between the /s and gs entering the other A$ = amphtudes Hiq, -^2,0, and -^4,0- The 
functions / contributing to if2,o ^ the other hadron-hehcity conserving amphtude (hehcity 
combination A$ = 0, A/ = +1/2, A^ = +1, Aj = +1/2) - are obtained from those entering 
H30 by interchange of the Mandelstam variables s ^->- u and the momentum fractions xi <-*■ 

yi- 



c{Q,D) _ AQ,D) 



fti:uH}-f{t^U-i} («"".-'"^.) 



ciq,D) _ Aq,D) 






,.(5) _ AS) 



fr 1 n = — fr 1 n (s ^^ M, Xi ^-> Ui) . 

•'{0,+i;+l,+i} ^{o,-i;+i,-i} ^ ^ V yv) 



(3.9a) 



(3.9b) 



(3.9c) 



The analytical expressions for the hadron-helicity-flip amplitudes -ffi,o, -f^4,0; -f^2,±i, and 
-f^3,±i are of similar length as the expressions quoted above and can be obtained from the 
authors on request. Like in Eqs. ( p.6[ ) the vector 3-point contributions to i/2,±i and -f^3,±i are 
proportional to the corresponding scalar ones. It is also interesting to observe that due to 
the (anti) symmetry properties of the Feynman graphs under interchange of the momentum 
fractions Z\ and Z2 (cf. Eqs. ( ^Tsj) ), only diagrams in which the photon couples to the s 
quark contribute to the A$ = amplitudes. 

Actually we have computed all the functions / and g, and relations like ( p.8|) and ( p.9|) 
served as a check of our analytical expressions. Furthermore, we have verified t/(l) gauge 
invariance with respect to the photon and S'[/(3) gauge invariance with respect to the gluon. 
The proof of gauge invariance is facilitated by the observation that not only the sum of all 
16 tree diagrams provides a gauge invariant expression, but rather each of the functions / 
and g is gauge invariant by itself. Finally, we have recalculated a few of the diagrams by 
hand to confirm the outcome of Mathematica. 
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IV. NUMERICAL RESULTS 

Fig. 1^ shows the diquark model predictions (sohd hne) for da/dt['yp -^ $p] at p^ab — 
6 GeV together with the outcome of the pomeron- exchange mechanism (dotted hne) |^ 
and resuhs of a two-gluon-exchange model (dash-dotted and dashed line) |0,0 • A direct 
comparison of our predictions with experiment is not possible yet since data are only avail- 
able at low momentum transfers (t ^ 1.5 GeV^) where the perturbative photoproduction 
mechanism is certainly not the dominant one. The experimental situation, however, will 
hopefully improve as soon as the data analysis of the CEBAF-93-031 experiment will be 
completed. Nevertheless, Fig. || shows that the diquark model provides a reasonable extrap- 
olation of the low-t data. Thereby one should keep in mind that the cross section which 
decays exponentially in forward direction is expected to flatten around 90° (t and u large). 
The few available data for photoproduction of pseudoscalar mesons and also photoproduction 
of p and u mesons which extend to larger values of t exhibit just such a behavior [^ . As can 
be seen from Fig. |^ the forward cross section for production is reasonably well reproduced 
by simple Pomeron phenomenology (dotted line). As a mechanism for photoproduction 
of vector mesons Donnachie and Landshoff ||^ proposed that the photon fluctuates into a 
quark-antiquark pair which subsequently recombines to a vector meson. The resulting quark 
loop is connected to a quark in the proton via the exchange of a single pomeron which cou- 
ples to quarks like an isoscalar photon. The QCD-inspired version of the Pomeron exchange 
of Donnachie and Landshoff P] , in which the Pomeron is replaced by two non-perturbative 
(abelian) gluons, has been applied by Laget and Mendez-Galain [^ to photoproduction of 



$ mesons (dash-dotted line). The latter authors consider the two-gluon-exchange model as 
a link which connects the diffractive with the hard-scattering domain. At t ^ 2.5 GeV^ the 
two-gluon-exchange cross section exhibits a characteristic node which can be understood as 
an interference effect between the two Feynman diagrams which enter the photoproduction 
amplitude. However, this picture changes drastically if the two gluons are allowed to couple 
to different quarks in the proton. Recently Laget |3^ tried to estimate such contributions 
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by modeling quark-quark correlations inside a proton via a simple wave function. He found 
that the node in the cross section is completely washed out and that the contributions in 
which the two gluons couple to different quarks inside the proton start to dominate the 
photoproduction of $ with increasing t (dashed line). Actually, for |t| ^ 4 GeV^ the result 
of Laget becomes comparable with the diquark- model result. This is not surprising since 
we know from the perturbative analysis that the hard-scattering mechanism, i.e. contribu- 
tions from diagrams without loops in which all the hadronic constituents are connected via 
gluon exchange, should become dominant for large values of t and u. Correlations between 
hadronic constituents are automatically accounted for within such an approach and hence 
also within our diquark model. 

To check whether our model provides a reasonable energy dependence we have also 
calculated the differential cross section at a much larger photon energy (pj^j^ = 2611 GeV). 
The result is depicted in Fig ^ (solid line) along with a two-parameter fit of |t| < 0.5 GeV^ 
data from ZEUS |3^. The hatched area indicates the uncertainties in the forward cross 
section and the slope parameter. Keeping again in mind that the cross section is expected 
to flatten with increasing t the magnitude of our prediction occurs still to be within the range 
of an extrapolation of the low-t data. The dashed line in Fig. ^ represents a leading-order 
calculation in which mass effects are neglected. This means that only the two hadron- 
helicity conserving amplitudes if2,o and H^q are taken into account in the leading-order 
cross section. A comparison between the full calculation (solid line) and the leading-order 
calculation (dashed line) reveals that mass effects are still sizable and amount to ~ 25-40% 
of the full differential cross section. They increase in forward and backward direction and 
are smallest around 90°. 

At this point we want to comment on another attempt to calculate photoproduction of $ 
mesons within the diquark model. The authors of Ref. ||2^ use a different parameterization 
of the model. Only the hadron-helicity-conserving amplitudes iJ2,o and H^^q are taken into 
account in their calculation. With their numerical method they encounter difficulties in 
calculating the real parts of these amplitudes. But they assert that the real parts are 
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small for their parameterization and therefore neglect them. Their results for differential 
cross sections are about one order of magnitude smaller than ours. This is plausible if one 
looks at their values of fg and fy, i.e. the "normalization" of the quark-diquark DAs, and 
the cutoff masses in the diquark form factors. Both are considerably smaller than ours. 
In contrast to the authors of Ref. PB[ we think that the present experimental situation 



in $ photoproduction does not allow to discriminate a particular parameterization of the 
diquark model. However, we want to point out that our parameterization has already 
been applied successfully to a variety of other exclusive processes [|l5H l9||. On the other 
hand, it is still unclear whether these processes can also be described satisfactorily with the 



parameterization of Carimalo et al. p6|. Their parameterization has been determined by 
means of the proton magnetic form factor. Apart from $ photoproduction it has only been 
applied to charmonium decays into the pp'y final state l^^ with the result that the data 
could not be reproduced. 

In the kinematic region, where t and u are large, the differential cross section at fixed 
scattering angle is expected to exhibit a scaling behavior of approximately s~^, which is 
supposed to indicate that the hard photoproduction mechanism becomes relevant. Fig. |^ 
shows s^ da/dt for photon energies of pj^^j, = 5, 10, and 20 GeV. 5 GeV is presently the 
upper limit of CEBAF, 10 GeV will be reached by the planned upgrade of CEBAF and 20 
GeV is a value which could be reached at future facilities like ELFE or a further extension 
of CEBAF which are presently still under discussion. Comparing the diquark-model results 
for these energies one observes that one is still far away from the s~'^ scaling behavior. The 



reasons have already been mentioned in Sec. III. A closer inspection shows that the cross 



section at 9cm = 90° scales rather like s~^ than s""^ for 10 < s < 40 GeV. Starting with 
point-like diquarks and a constant ag one would expect that the leading-order cross section 
(without mass effects) should scale like s~^. By including the diquark form factors and cross- 
section contributions from hadron-helicity-fiip amplitudes, which are already suppressed by 
one power of s as compared to the leading-order contributions, the scaling behavior changes 
already to nearly s~^. Finally, the approximate s~^ decay is caused by the fact that Og is 
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taken as a running coupling constant (with argument (iu/s)). 

Fig. P shows how the differential cross section (at p/ab — 6 GeV) is inffuenced by different 
contributions within our model. A comparison of the solid (full calculation) and the short- 
dashed (mass effects neglected) line reveals, e.g., the influence of mass effects. For pj^j^ = 
6 GeV and cos(^cm) = 90° they amount to ~ 40% of the full cross section. The dash-dotted 
line represents the pure S-diquark part of the cross section. It varies between 15% and 30% 
if one goes from backward to forward angles. Since S- and V-diquark contributions add 
coherently (and constructively, cf. Eq. ( |3.6| )) this means that the corresponding amplitudes 
are of the same order of magnitude. The differential cross section for the photoproduction 
of longitudinally polarized $s is given by the long-dashed curve. It amounts to ~ 60% 
of the full cross section at cos(^cm) = 90°. This fraction increases in backward direction 
and decreases in forward direction. Finally, we have also examined the influence of the 
exponential in the $ DA (cf. Eq. (|2.10| )). Neglecting this exponential one ends up with the 
dotted curve which is nearly indistinguishable from the full calculation. This also shows that 
we are not very sensitive to the 2:1 ^ 0, 1 endpoint regions when performing the convolution 
integral Eq. (|2l| ). 

As examples for single-polarization observables the beam asymmetry S^. and the target 
asymmetry Ty are depicted in Fig. |^. The beam asymmetry would be already non-zero, 
if only the leading-order helicity amplitudes if2,o and H30 were taken into account (cf. 
Eq. ( p. 3D ). It would not even be necessary that the helicity amplitudes acquire a non-trivial 
phase. On the other hand, a non-trivial target asymmetry demands for both, non-vanishing 
hadron-helicity-fiip amplitudes and non-trivial phases of the amplitudes (cf. Eq. (|3.4])). As 
already mentioned at the end of Sec. [11 B| phases are generated through propagator singu- 
larities, hadron-helicity-fiip amplitudes occur in our model as mass corrections. As it turns 
out within our model, also the beam asymmetry is considerably affected by the inclusion of 
the hadron-helicity-fiip amplitudes. In a recent paper [Q the influence of nucleonic reso- 
nance effects on polarization observables in $-meson photoproduction has been investigated 
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within a constituent-quark model in which the diffractive contribution is produced by t- 
channel Pomeron exchange. It is interesting to observe that the resuh of these authors for 
the target asymmetry at a rather small photon energy of pj^j^ = 2 GeV is comparable in size 
to our result at p^^ — 6 GeV. This is remarkable since two completely different concepts 
are at work: on one hand nucleonic resonance contributions, on the other hand a purely 
perturbative mechanism. Also for the beam asymmetry some similarities can be found. The 
corresponding predictions of both approaches are again comparable around 6*011 = 90°. A 
(somewhat smaller) negative beam asymmetry has also been found in another model which 
is based on pomeron exchange, pion exchange, and a direct knockout production from the 



strange-quark sea in the proton [^. It would be interesting to see experimentally, whether 
polarization observables like the target asymmetry Ty, the recoil polarization Vy', or the 
vector-meson polarization Vy', which vanish in leading-order perturbative QCD, are still siz- 
able at momentum transfers of a few GeV. Such observables would be a more sensitive tool 
than the unpolarized differential cross section to figure out whether the kind of nonpertur- 
bative ingredients (diquarks, effective constituent masses) which are taken into account in 
our perturbative approach, are already sufficient, or whether much more non-perturbative 
physics is still at work at a few GeV of momentum transfer. 

At the end of this section we want to comment on the perspective to extend this calcu- 
lation to photoproduction of J/^. It is a simple matter, since we have already assumed a 
flavor dependent DA for the vector meson (cf. Eq. ( p.lOp ). The only thing to do is to replace 
the $ mass by the J/\E' mass and the decay constant /* by Z"^/*. One, however, has to keep 
in mind that we have applied an expansion in the parameter (rriH/yS). In order that this 
parameter is approximately of the same size in $ and J/\l/ production the photon lab energy 
has to be increased approximately by a factor {mji^i/inqyY = 9-2 when going from $ to J/\&. 
A second criterion for the reliability of our approximate treatment of mass effects is that the 
coefficient of the first mass correction term should become small compared to the leading 
order. This coefficient is angular dependent and becomes large at small scattering angle. 
That means, if the size of the mass corrections relative to the leading order should be about 
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the same in $ and J/^ production {{mH/ys) comparable), one has to consider much larger 
momentum transfers for J/\l' than for $ production. A more detailed discussion of this 



point along with some numerical results can be found in Ref. |^0[- ^ possible way to apply 
the perturbative approach to J/\E' production at a few GeV of momentum transfer would be 
to retain our prescription for the choice of the effective masses, but to refrain from the final 
expansion in (m/^/Vs). This has not been attempted yet since the analytical expressions 
for the amplitudes become very lengthy and the treatment of propagator singularities in the 
convolution integral, Eq. ( p.l|) , becomes much more involved. 

V. SUMMARY AND CONCLUSIONS 

We have investigated the reaction 7p -^ $p in the few-GeV momentum transfer re- 
gion where unpolarized differential cross-section data from the CEBAF-93-031 experiment 
at JLAB are expected to become available very soon. Our theoretical approach is based on 
perturbative QCD supplemented by the assumption that baryons can be treated as quark- 
diquark systems. The same approach has already been successfully applied to other photon- 
induced hadronic reactions |T5|-|19|. The motivation for introducing diquarks is, above all, to 
extend the perturbative approach from large down to moderately large momentum transfers. 
There are mainly two reasons why this can be achieved with diquarks. First, the overall 
momentum transfer has to be shared between less hadronic constituents than in the pure 
quark HSA. Hence the gluons which keep the hadronic constituents collinear are on the 
average harder and thus the running strong coupling constant becomes smaller. The second 
reason is that a certain class of non-perturbative effects, namely quark-quark correlations 
inside a baryon, are effectively accounted for by means of diquarks. The fact that a heavy 
quarkonium is produced in the final state entails a considerable reduction in computational 
effort. Only a small fraction of the Feynman diagrams which contribute in general to photo- 
production of arbitrary mesons occurs in $ production. To be more specific, these are just 
those diagrams in which the photon fluctuates into the heavy quark- antiquark pair which is 
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then connected to the quark and the diquark in the proton via the exchange of two gluons, 
respectively. Both, the scalar and the vector-diquark component of the proton have been 
taken into account. We have also improved the diquark model in the sense that hadron-mass 
effects in the hard-scattering amplitude have been included in a systematic way. This has 
been accomplished by assigning to each hadronic constituent an effective mass xtuh {x . . . 
longitudinal momentum fraction of the constituent) and expanding the scattering ampli- 
tude in the small parameter [mH/ys). Only the leading and next-to-leading order terms 
in this expansion have been kept. The leading-order term provides for the hadron-helicity- 
conserving amplitudes and agrees, of course, with the massless result. The next-to-leading 
order term represents a power- correction and gives rise to non- vanishing hadron-helicity-flip 
amplitudes. It enables us to make predictions for polarization observables which require the 
flip of a hadronic helicity like, e.g., the target asymmetry or the recoil polarization. We want 
to point out that our treatment of mass effects still preserves gauge invariance with respect 
to the photon and the gluon. 

Our numerical studies have been performed with the diquark-model parameters and the 
quark- diquark DAs proposed in Ref. JT^. The DA of the $ meson has been taken from 
the literature |2^. It fulfills QCD sum- rule constraints. Hence we have been able to make 



predictions without introducing new parameters. We have paid special attention to the 
correct and numerically robust treatment of propagator singularities. In the absence of 
large-t data (t > 4 GeV^) our results can only be compared with the available data which 
go up to t ~ 1.5 GeV^. Our predictions for the differential cross section at two rather 
different photon energies (pjl^j, = 6 GeV and 2600 GeV) which are typical for CEBAF and 
HERA, respectively, occur to be still within the range of an extrapolation of the existing 
low-t data. It is worth mentioning that the predictions of a two-gluon-exchange model, in 
which the gluons are considered as a remnant of the pomeron [^, become comparable to 
our result for t ^ 4 GeV^ if the gluons are allowed to couple to different quarks in the 
proton. This is not surprising because contributions in which all the hadronic constituents 
are (approximately) kept coUinear by means of gluon exchange to prevent the hadrons from 
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breaking up become increasingly important as soon as one enters the hard-scattering domain. 
Just those contributions are assumed to be the dominant ones within our HSA-based diquark 
model. We have also investigated within our model whether the (asymptotic) s~^ scaling 
behavior of the differential cross section at fixed angle is already realized for s values between 
10 and 40 GeV. After dividing out the strong running coupling (which occurs to the sixth 
power in the cross section) an approximate s~^ scaling has been found. Our investigation of 
hadron-mass effects has revealed that they are still non-negligible in the few- GeV momentum 
transfer range and photon energies of the order of 10 GeV. Correspondingly, also polarization 
observables which require the flip of a hadronic helicity, like the target polarization, occur 
to be sizable. Finally, we have considered the possibility to extend the present calculation to 
photoproduction of J/^ mesons. In order that our treatment of mass effects remains reliable 
for J/\E's one has to go to much larger energies and momentum transfers where experimental 
data cannot be expected. The situation with respect to the momentum transfer could 
perhaps be improved by refraining from the final expansion in (m///vl)- This, however, will 
considerably complicate the treatment of propagator singularities and make the analytical 
expressions for the amplitudes much more complicated. Therefore it has not been attempted 
as yet. 

Our diquark-model calculations for '-/p -^ $p will meet the first severe test in the near 
future when the data analysis of the CEBAF 93-031 experiment will be completed. With 
these data, which go up to t values of about 5 GeV^, one will get a first glimpse how 
different production mechanisms (pomeron exchange — ^ gluon exchange — ^ hard-scattering 
mechanism) in high-energy photoproduction evolve with increasing momentum transfer. 
However, in order to really touch the hard-scattering domain much higher photon energies 
than those which can be reached at the present stage of CEBAF would be desirable. In 
addition to unpolarized differential cross sections, polarization observables which require 
the flip of hadronic helicities, could serve as an indicator for the relative importance of 
hard and soft photoproduction mechanisms, since hard scattering is closely connected with 
hadronic helicity conservation. 
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FIGURES 
FIG. 1. Graphical representation of the hard-scattering formula, Eq. ( |2.1| ), for 7p — > M p. Xi, 

Hi, and Zi denote longitudinal momentum fractions of the constituents. 

FIG. 2. Eight of the, altogether, sixteen tree-level diagrams which contribute to the elementary 
process 'jqD -^ QQqD. The remaining diagrams are obtained by interchanging the two gluons. 
We have also indicated how the diagrams are grouped to gauge invariant expressions. 

FIG. 3. Differential cross section for jp -^ ^p versus t at p^j^ = 6 GeV. The solid line 
corresponds to the diquark-model prediction. Cross sections resulting from the Pomeron-exchange 
mechanism [Q| (dotted line) and a two-gluon-exchange model |3l| (dash-dotted line) are plotted for 
comparison. The dashed line represents an attempt to extend the two-gluon-exchange model by 
including contributions in which the two gluons couple to different quarks inside the proton |32|. 
Data are taken from Ballam et al. [33| and Barber et al. |34]. 



FIG. 4. Differential cross section for jp ^ ^p versus t at p^^^ = 2611 GeV {s = 4900 GeV^). 
The diquark-model prediction with (solid line) and without (dashed line) mass effects is compared 



with a parameterization of the low-t ZEUS- data |36]. The hatched area indicates the uncertainties 



in the forward cross section and the slope parameter. 

FIG. 5. Differential cross section for jp -^ ^p scaled by s'^ versus cos(^cm)- Diquark-model 
results for p'^^-^ = 5 (solid line), 10 (long-dashed line), and 20 GeV (dashed line). 

FIG. 6. Differential cross section for 7p — > ^p scaled by s^ versus cos(6'cm) at p^^^ = 6 GeV. A 
comparison of different contributions within the diquark model: full calculation (solid line), mass 
effects neglected (short-dashed line), S diquarks only (dash-dotted line), longitudinally polarized 
$s only (long-dashed line), exponential in <I>-DA neglected (dotted line). 

FIG. 7. Diquark-model predictions for beam (S^.) and target (Ty) asymmetry in ^p — > ^p 
scattering at p^^-^ = 6 GeV. 
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C.F. Berger and W. Schweiger: Hard exclusive photoproduction of O mesons, Fig.l 
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C.F. Berger and W. Schweiger: Hard exclusive pliotoproduction of <I> mesons, Fig.2 
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C.F. Berger and W. Schweiger: Hard exclusive photoproduction of <I> mesons, Fig. 4 
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